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We investigate the inelastic spin-flip rate for electrons in a quantum dot due to their contact 
hyperfine interaction with lattice nuclei. In contrast to other works, we obtain a spin-phonon 
coupling term from this interaction by taking directly into account the motion of nuclei in the 
vibrating lattice. In the calculation of the transition rate the interference of first and second orders 
of perturbation theory turns out to be essential. It leads to a suppression of relaxation at long 
phonon wavelengths, when the confining potential moves together with the nuclei embedded in 
the lattice. At higher frequencies (or for a fixed confining potential), the zero-temperature rate 
is proportional to the frequency of the emitted phonon. We address both the transition between 
Zeeman sublevels of a single electron ground state as well as the triplet-singlet transition, and we 
provide numerical estimates for realistic system parameters. The mechanism turns out to be less 
efficient than electron-nuclei spin relaxation involving piezoelectric electron-phonon coupling in a 
GaAs quantum dot. 
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I. INTRODUCTION 

Electron spin dynamics in mesoscopic devices has been 
attracting a lot of attention recently in the context of 
spintronicsi and quantum computations*^. A crucial fea- 
ture of this dynamics is the relaxation of the electron's 
spin due to the interaction with an environment. Gen- 
erally speaking, the coherence of an electronic spin state 
vanishes during the time T2, which limits the possibil- 
ity of coherent manipulation of qubits, while relaxation 
to thermal equilibrium occurs during another time Ti, 
which is usually larger than T-^. 

Several types of spin-dependent interactions can give 
rise to electron spin relaxation, e.g. the electron spin- 
orbit interaction^i'^i'^i^iQi^°i"i^'^i^^i^^ and the electron- 
nuclei hyperfine interaction'^i^^i^^i^'^i-^^i^^i^°i^^. Their ac- 
tion depends essentially on the dimension of the system. 
Systems of zero dimension, i.e. quantum dots (QDs), 
are characterized by a discrete electron energy spectrum. 
In this case energy conservation in the spin-flip process 
usually can be fulfilled only by transferring the energy to 
another subsystem, e.g. phonons. The energy transfer 
includes both the Zeeman energy of the electron spin in 
an external field and possibly the energy of an orbital 
transition. A discussion of other electron spin relaxation 
mechanisms not mentioned above and relevant to QDs 
can be found in Ref. ITol for example. Electron spin re- 
laxation in a QD due to hyperfine interaction alone, in 
the absence of an ext ernal m agnetic field, was investi- 
gated recently in Refs. I22l23i 

Many years ago, the phonon-assisted electron spin-fiip 
transition between Zeeman sublevels due to hyperfine 
interaction with an impurity nucleus and lattice nuclei 



was considered for impurity-bound electrons in silicon^Si, 
where the authors investigated the nuclear polarization 
in the Overhauser effect. The process that has been con- 
sidered in that paper is associated with a crystal dilation 
and a corresponding adiabatic change in the electron ef- 
fective mass and, as a consequence, in the electron enve- 
lope wave function and the hyperfine coupling constant. 
Recently, electron spin relaxation due to the hyperfine 
contact interaction has been readdressed with an appli- 
cation to GaAs QDsiSiiS. The transition amplitude was 
calculated in second-order perturbation theory, describ- 
ing the action of the hyperfine contact interaction and a 
spin-independent piezoelectric electron-phonon coupling. 

In this paper we will analyze another spin relaxation 
mechanism provided by the combination of hyperfine 
contact interaction and the influence of phonons on the 
electron inside a QD. In our approach we take into ac- 
count directly the phonon-induced motion of nuclei which 
are coupled to the electron spin through the hyperflne in- 
teraction. The electron-phonon interaction appears via 
the displacement field shifting the positions of the nuclei, 
and therefore is independent of the piezoelectric coupling 
that applies only to crystals without inversion symmetry 
(such as GaAs, but not Si). Moreover, this mechanism al- 
lows the electron-nucleus spin fiip-flop with a simultane- 
ous emission of a phonon to appear already in first-order 
perturbation theory. Nevertheless, it will turn out that 
it is necessary to keep as well the second-order terms as- 
sociated with the motion of the electron confining poten- 
tial, since they lead to a crucial cancellation of first-order 
terms at low frequencies of the emitted phonon, thereby 
suppressing the relaxation rate. The physical reason be- 
hind this is the following: Long-wavelength phonons dis- 
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place both the lattice nuclei and the electron's confining 
potential in the same way. However, it is only the rela- 
tive motion of the electron with respect to the nuclei that 
enables a transition. Therefore, the influence of long- 
wavelength phonons is suppressed. If, on the other hand, 
the confining potential can be considered as fixed or it 
moves independently from the lattice nuclei, this sup- 
pression does not apply any more, and the destructive 
interference between first-order and second-order terms 
is broken. 

The present article is organized as follows: In Section 
UTI we introduce the Hamiltonian of our model, includ- 
ing the effect of lattice vibrations on the hyperfine cou- 
pling and the confining potential. After that, we derive 
the electron spin fiip transition amplitudes due to this 
perturbation, discussing the partial cancellation of terms 
and the dominating contribution. We calculate the tran- 
sition rate, Eq. (|23|l . for the case of Zeeman-split sub- 
levels of the electron ground state (Sec. IIII|I . and perform 
an analogous derivation for the case of a triplet-singlet 
transition (Sec. II V|) . Finally, in Section El we look at 
numerical estimates for realistic system parameters. 



II. THE MODEL 

In our model, we will assume the displacement of the 
QD confining potential to be described by the phonon 
displacement field evaluated at the center of the dot 
(which we take to be the origin). Note that this is anal- 
ogous to the case of an impurity-bound electron. Any 
more detailed description (e.g. allowing for a distortion 
of the potential) would require further specifications con- 
cerning the way this potential is applied to the dot, but 
would not add significantly to the realism of the present 
model. 

In the effective mass approximation the Hamiltonian 
of the system of electrons and phonons in the QD has 
the following form, if the perturbations due to electron- 
phonon and hyperfine interactions are excluded: 



P7 + V{r^ 

2m* 



g*fj,BSi ■ B 



,(1) 



where = — i?iVri +eA-(ri)/ c is the kinematical momen- 
tum operator. This Hamiltonian describes interacting 
electrons with effective mass m* and effective g-factor g* , 
localized in a static external potential V{y), in the pres- 
ence of a magnetic field B (with a corresponding vector 
potential A(r)), and free phonons. We neglect the Zee- 
man splitting of nuclear spin states, given the small value 
of the nuclear magnetic moment. 

The hyperfine contact interaction of electrons of spin 
Si at positions with nuclei of spin Ij at positions Rj 



has the form [J]: 

j 

= Av^Y.^,-i,5{v,-K,), (2) 

j 

where the hyperfine coupling constant is determined as 



A 



31 vq 



■V, 



with 77 = |u(0)P being the square of the Bloch amplitude 
at the site of the nucleusi^, /is > the Bohr magneton, 
fii the nuclear magnetic moment, and vq the size of the 
unit cell. The factor 77, which is usually of the order of 
10^ to 10*^, depending on the material, makes the Fermi 
contact term much more efficient than the other terms of 
the electron-nucleus hyperfine interaction-. For a more 
complicated lattice structure the hyperfine constant A is 
defined as a sum of the individual constants Aj over all 
the nuclei in the unit cell, taking into account different 
values of spins and parameters 77 for different nuclei^^. 

The positions of the nuclei deviate slightly from equi- 
librium, due to the lattice vibrations: Rj = R^ -I- u(Rp. 
The lattice displacement field u(R°) is described via the 
phonon creation-annihilation operatorsSi, 



k,A 



{bk,\ + b\ ^) ek,A e 



ik-R" 



(4) 



where ek,A is the polarization vector of a phonon with 
wave vector k in branch A, uJk.\ is its frequency, Vph is 
the volume of the crystal in which phonon modes are 
quantized, and p is the crystal mass density. 

As a consequence, the total Hamiltonian acquires the 
following term which can lead to a nucleus-electron spin 
flip-flop process combined with the emission of a phonon: 

VHF-ph{r^) = ~AvoJ2 • h • VrJ S{r, - RO). 

j 

(5) 

In the following, we will omit the index for nuclear 
equilibrium positions. 

Likewise, the vibrations of the confining potential are 
described by 

E'5nr.) = -E("(0)-^rjnr.) 

i i 

= u(0) [i/o,E ^P'] = -^"(0)E t^O' t^O'''^]] (6) 

i i 

We note that the total electron momentum commutes 
with the electron-electron interaction potential. 

Thus, the total perturbation to the Hamiltonian Hq^ 
Eq. (Q, is given by the three terms described above: 

<5i? = E [^^^""i) + yHF{vi) + VHF-phiri)] . (7) 
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III. TRANSITIONS BETWEEN ZEEMAN 
SUBLEVELS 

At first, we will consider the transition of a single 
electron between the Zeeman-split spin levels of the QD 
ground state. 

The initial state of the system is given by the direct 
product of electron, nuclear and phonon states: 



|«7v) \iph)- 



(8) 



Let n denote the direction of magnetic field. The initial 
electron state \ie) = ® |n_) is given by the product 

of the spin wave function |n_), which is an eigenfunction 
of the equation (n • S)|n±) = ±i|n±), and the properly 
normalized ground state coordinate wave function 



1 



1 
V 



(9) 



Here V is the effective volume of the dot. In GaAs, where 
the electron g-factor is negative, the state |n_) corre- 
sponds to the maximum of energy. The initial nuclear 
spin state \iN) is a direct product of states of all indi- 
vidual nuclei. We will average over the initial phonon 
field state in the end, by inserting mean phonon occu- 
pation numbers given by the Bose distribution function 
n„ = l/(e'*"/fci'^- 1). 
In the final state 



I/) = l./e) ® \fN) ® \fph), 



(10) 



the electronic spin points into the opposite direction, 
l/e) = IV'(r)) ® |n+). Nuclear, I/at), and phonon, \fph), 
final states are determined by the action of the pertur- 
bation potential VuF-ph which changes the electron and 
nuclear spin states while conserving the total spin of the 
electron-nuclei system (flip-flop process) and creates a 
phonon with energy corresponding to the energy differ- 
ence between electron initial and final states. 

The corresponding transition matrix element is pro- 
vided by first order perturbation theory in the potential 
VHF-ph{r)- 

(/|^i?|z)W 

^-Y.^fph\n{'^j)\iph){feM'^rVhs{v-B.j)\l,^M)-{ll) 



where the sum is over all intermediate states which dif- 
fer from the initial and final states, and iJ^"^) refers to 
electron energies only. In writing down the energy de- 
nominators we have used the fact that the hyperfine per- 
turbation only changes the electronic energies, and that 
initial and final total energies will be the same. Note that 
the complete electronic energy includes the Zeeman en- 
ergy as well, and the difference between initial and final 
electron energies is accounted for by the energy of the 



emitted phonon, Tiuj = e[^^ — E 



pie) 
f 



Now we employ the expressions for the perturbation 
SV{r), Eq. IHI, in order to rewrite the previous formula 
in a form that displays the relation to the first order 
amplitude: 

(/|,5i?|i)(2) = (fp,,\uiO)\iph) [{fe,N\^rVHFir)\le,N) 
{fe,N\ r \me,N) {me.N\VHFi'r)\ie.N) 



* 2 

m uj 



(/e.jv|VH"F(r)|me,Af)(me,Af| r \ie,N) 
pie) _ pie) 



(13) 



We have used the expressions for 5V{r) in terms of the 
commutator of Hq with both the kinematical momentum 
and the position operator, as well as the fact that Vhf 
commutes with the position. 

The total amplitude of the transition is the sum of the 
terms Ijllf) and H13|l . We will regroup it into two parts. 
The first one consists of Eq. (|11() and a contribution of a 
similar form, the first term of Eq. (|13|l . 



^1 - E (fphim - "(R^- )M ife-NlVrVkf (r-R, ) I 
3 



le,N) 



(14) 

It contains the difference between the phonon displace- 
ment fields evaluated at the origin and at the nucleus 
position, respectively, which is analogous to taking the 
divergence of the displacement field. This difference van- 
ishes as I kRj I for long- wavelength phonons and describes 
the important cancellation which is only found if first and 
second orders of perturbation theory are combined prop- 
erly. 

The second part. 



It is important that an alternative process is possible 
for the transition between the^ same two states, where 
the remaining two terms in 6H (namely 5V and Vhf) 
contribute in second order perturbation theory, yielding 
an amplitude that is of the same order of magnitude as 
Eq. (ini): 



{f\6m'^ = E' 



{f\6V{r)\m){m\VHF{r)\i) 



{f\VHF{v)\m){m\5V{v)\i) 



E\ 



ie) 



(12) 



A2 = -{fphHmph) 

(/e.A'l r \me^N){meMVHF['r)\ie,N) 



X m*Lj^ 



E, 



ie) 



E: 



ie) 



{fe,N\VHF{r)\me,N){me,N\r \ie,N) 



E 



ie) 
f 



E: 



ie) 



(15) 



contains a sum over intermediate states. For a single elec- 
tron making a transition between the Zeeman sublevels 
of its ground state, we can set Ef^ —eIji' = £Q — i£m — fi^) 
in the first sunimand and — Em^ — (sq — Jiuj) — Em in 
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the second. Note that in the two sums in Eq. (|15(l inter- 
mediate states with the same orbital energies differ 
by the Zeeman energy. 

In order to render the following discussion concrete, 
we will now specify the confining potential explicitly. 
We consider a QD which is formed in a two-dimensional 
electron gas (2DEG) by an external symmetric parabolic 
potential. The confining potential in z-direction is usu- 
ally modelled by a square well in vertical QDs and by 
a triangular-shaped potential in lateral QDs. We ne- 
glect the contributions from higher excited states in the 
z-potential, given their large energetic separation, and 
restrict our discussion to the ground state Xo{z)- In the 
presence of an external magnetic field perpendicular to 
the X — y-plane, the electron wave functions in the lat- 
eral dimension become the Darwin-Fock solutions (see 
e.g. Ref. 125^ . with the effective confining frequency 



ujQ = y^n^^+ujf/A, where ilo is the strength of the 
parabolic potential and ojc — eB±/{m*c) the cyclotron 
frequency in an external magnetic field with the compo- 
nent B± perpendicular to the layer. The energy spectrum 
for these states is e„A — (2»t. -I- |A| -t- l)hujo — Xhujc/2. 

In a harmonic oscillator the coordinate vector induces 
transitions only between nearest orbital levels. Thus, the 
sum in Eq. IjlSII reduces to two terms only correspond- 
ing to the transition from the ground state (j)oo to states 
(j)nX, with n = and A = ±1 being the radial and the 
angular momentum quantum numbers, respectively. For 
example, in cylindrical coordinates we rewrite the scalar 
product of the phonon polarization vector and the coor- 
dinate vector as er = (e+pe^"^ + e-pe^'^)/^/2 and note 
that (?!)ooPe^"^/\/2 ~ l(t>o±i, where I — ^Jfi/{mFujo) is 
the length scale which determines the spatial extent of 
the electron wave function in the parabolic well in the 
presence of an external perpendicular magnetic field. 

We now compare the two parts of the total amplitude 
that are given by Eqs. (|14|l and (|15() . First, we note that 
the matrix element (/e.Ar|VrV/i/(r — Rj)|ie.Ar) in the ex- 
pression for Ai is proportional to |VR,(/)o(Rj)|, where the 
gradient in the z-direction can be estimated as ~ 1/2:0, 
with zq the transverse dimension of the dot. Gradients 
in the lateral directions are smaller by a factor zq/1 ^ 1. 
The Zeeman splitting, Ulj — \g* ^bB\, is much less than 
the orbital energy splitting, Tiuj <^ hojo, in typical experi- 
ments on spin relaxation in QDs. Thus, up to a common 
prefactor, we can use the following estimates for the ex- 
pressions of Eqs. (|14|l and (|15|l : 



lAiloc 



i[kl,l] 



A/s,l] 



, .2/ ' 

(16) 

where s is the mean sound velocity. Finally, we find that 
the ratio 1^2/^1 1 is suppressed by a factor whose form 
depends on the emitted phonon wave-length: 



Zq to S 

I loq loqI 



< 1 for fc/ < 1 



) «1 for fcZ»l. (17) 
I \ujoJ 

Taking this into account we will neglect A2 in what fol- 
lows. 

Now we rewrite the total transition amplitude in the 
following, more explicit form (retaining only the main 
contribution) : 



(/|<5i7|i)(*°*°') 



A 



Vq_ 
V 



J2 (e* ""^ 



X ie^.x-Vn,)^^-^,), (18) 

where S"_ = (n-|_|S'"|n_) (a sum over a in Eq. (|18|l is 
assumed) . 

By means of Fermi's golden rule we obtain the follow- 
ing expression for the transition rate (including a sum 
over final states and a proper average over initial states): 



2n A^ 
W = 



, + l)^i^(R„R,OG. 



VR,.(^g(Rj)-VR^,^2(J^^.,) 



(19) 



where N = V/vq is the number of unit cells inside the 
dot volume V . In H19(l we have employed the sound wave 
dispersion law in the form Wk = ks, i.e. we have neglected 
the difference in the transverse and longitudinal sound 
velocities in summation over phonon polarizations. This 
simplifies our formulas but should not change appreciably 
our numerical results. 

For clarity, in Eq. (|19|) we have combined some expo- 
nential factors from Eq. H18|l into the following expres- 
sion: 



^^(R„R,0=4e"< 



k R, k R,v , ^ 
■ sin ^ sin ^ . (20) 



We have also separated all spin-related factors into the 
correlation function (cf. Refs. 18 19): 



Gj,. ^S%Sl_{lN\ifi'^,\lN)av, 



(21) 



where the subscript, aw, indicates averaging over initial 
nuclear spin states (in our case over a thermal distribu- 
tion). 

At temperatures much higher than ~ 10~^ K ^ which 
is the order of magnitude of the nuclear spin-spin inter- 
action, the nuclear spins are not correlated, i.e. Gjji — 
GSjj'. We suppose that there are no other sources of 
average nuclear polarization either. This means, in turn, 
that the interference terms in Eq. (|19|) . stemming from 
different nuclei, j 7^ j', vanish. 
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With the help of the usual commutation rules for spin 
components and the equality (n_|S|n_) — — ^n, we ob- 
tain for each nuclear spin the formula 



r - (i • nf + (i • n) 



(22) 



which results in a correlation function G = + l): 

provided the average nuclear spin is zero. 

Taking all of this into account, the following expression 
describes our main result for the rate of the electron-spin 
relaxation between Zeeman sublevels of the ground state, 
due to the hyperfine-phonon mechanism considered here: 



1 



1 



LUl{I+l) 



sin(fc|R'| 



fc|R'| 



(23) 



The overbar in Eq. (|23|l indicates an average over the 
positions R' of all nuclei in the dot. If the electron enve- 
lope wave function changes little on the scale of distance 
between the nuclei, then this is just a spatial average over 
an electron localization volume. In deriving Eq. (|23|l . we 
have used the identity 



dn^ sin^(k- ^) = 27r 



sin(2fc|^|) 



(24) 



We note that the relaxation rate for an individual nu- 
clear spin is obtained by dividing Eq. H23|) by N . 

IV. TRIPLET-SINGLET TRANSITION 

We proceed in the same way in order to calculate the 
transition rate for two electrons that initially reside in the 
lowest-lying triplet state and decay towards the ground- 
state singlet. We suppose that the Zeeman splitting (pro- 
duced by external and/or nuclear magnetic fields) can be 
neglected as compared to the orbital energy spacing that 
defines the transition energy for this process. 

The wave function corresponding to the initial spin- 
triplet electron state is 



|V'^'-(ri,r2))® iTr), 



(25) 



where the coordinate wave function is assumed to be 
given by a Slater determinant (i.e. neglecting correla- 
tions): 



(po(ri)i^i(r2) - v?o(r2)(/?i(ri) 
VV2 



(26) 



Here (pi (r) is the wave function of the first excited sin- 
gle electron orbital state normalized according to Eq. 
For concreteness, we choose it to correspond to the quan- 
tum numbers n — and A = 1 (in an external magnetic 
field the energy of this state is less than £0.-1)7 and we 



note that there is no term in the total Hamiltonian which 
couples directly states with A = ±1. We write for the spin 
part (as in Ref. 0) 



11, +1) 



V2 ' ' V2 

(27) 

where the coefficients i'x,y,z determine the initial su- 
perposition of degenerate states \S^m) with different z- 
components m = ±1,0 of the total spin 5 = 1 of two 
electrons. 

For the final spin-singlet state we have 



Si\ 



\4'^\v^,V2))®\Sl), 



(28) 



where \Si) denotes the singlet spin state (5* = 0), and 
the coordinate wave function is given by the expression: 



yo(i'i)¥^o(r2) 
V 



(29) 



Again, the transition amplitude is given by Eqs. 
and (|15|l . where we have to introduce sums over electron 

coordinates: e.g. VrT4/(r-R-j) ^ ^rit4/(ri -Rj), 
etc. 

In the present case, the energy of the emitted phonon, 
TiLj^ is equal to the single-particle energy splitting, TIliJq. 
In the expression for A2^ Eq. (|15|l . the contribution from 
the spin-singlet intermediate state dominates, due to the 
small denominator (18| given by the exchange splitting 
5 ST- However, although 5st is smaller than Tiljq^ it 
still has the same order of magnitude^^, and the ratio 
~ {zq/1)[ijJo/5st) is still much less than unity 
(here we have kl 3> 1). Therefore we can once more 
neglect A2. 

We obtain for the triplet-singlet transition rate: 



W 



ST ^ 



h 
1 r 



(5(?i^^k - hu;) 
, + l)^^^(R„R,OG|,T 



VR^.((^o(Rj)<y5i(R,)) • Vr^,((^o(R,')^i(R/)) 

(30) 

where the correlation function is: 



G 



ST 
if 



l^*''l^^{^N\It■I';\^N)aJJ 



-/(/+!)%. .(31) 



Again, the nuclear spin state is averaged over a com- 
pletely unpolarized thermal distribution. 
Finally, the relaxation rate in this case is: 



1 



1 



Tf^ N 6TThps' 



u;I{I + l) 



1 - 



sin(fc|R'| 



fclR' 



[Vr.(^o(R')^i(R'))] ■ (32) 
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V. NUMERICAL ESTIMATES AND 
DISCUSSION 

In order to estimate the rates in both of the cases that 
have been considered above, we take into account realis- 
tic dimensions of typical quantum dots. Usually the lat- 
eral length of the dot is much larger than its transverse 
dimension, I ^ zq, and we use the following approxima- 

tion: [Vr-(^o(R')^i(R'))]' ^ [^n'V^UH'f ^ 
(which we suppose holds for the average over nuclear po- 
sitions R'). For QDs with a disc shape, in the typical 
limit kzQ < 1, we can obtain a simple analytical expres- 
sion for the average, 



f{kl) 



sin(fc|R'|) 



A:|R'| 



1 



2(1 - cos(A:/)) 
(fc02 ' 



(33) 



which can be well approximated by min{(fc^)^/12, 1}. Wc 
remark that this factor, which is present in the case of 
the co-moving confining potential, is absent when the 
confining potential position is fixed. This is the only 
difference in the results for the relaxation rates in these 
two cases, and it becomes important only in the limit of 
small phonon frequencies, kl <^ 1. 

We can write the relaxation rate between Zeeman sub- 
levels and triplet and singlet electron states using the 
same approximate expression (the difference lies in the 
energy scale lu): 



Ti~ ?, N irfips^zl 



I{I+l)f{kl). (34) 



The linear dependence of the relaxation rate on the 
phonon frequency (for phonon wave lengths smaller than 
the size of the QD, i.e. for sufficiently strong magnetic 
fields in the case of Zeeman sublevels relaxation) sets 
our mechanism apart from those considered earlieriiSiiS, 
where a cubic phonon frequency dependence is expected 
for low temperatures, ksT <ti huj (and a quadratic one 
in the opposite limit). In addition, our result does not 
depend either on the proximity of the nearest level, in 
contrast to the spin relaxation rates calculated in Refs. 

mia 

The result for the electron-nuclei flip-flop transition 
rate obtained in Ref. [3 is larger than our result, Eq. 
(|34|l . by a factor ~ (jzo/l)'^ (provided kl ^ 1), which 
is of the order of 10 (7 ~ 50) for typical QD sizes cited 
below, but can be smaller for large QDs. We should be 
cautious, however, in applying directly the reasoning of 
Ref. TS to the triplet-singlet transition, when the emit- 
ted phonon energy corresponds to the electron binding 
energy and hence the condition of adiabatic electron mo- 
tion in a vibrating lattice (used in that work) will not be 
fulfilled. 

In GaAs all nuclei have spin I — ^, and A is the sum 



of Ai over all the nuclei in the unit cell: A 



~ 90 

peV-i^. The mass density is p ~ 5.32 x 10^ kg/m^, and 
we approximate the mean sound velocity by the velocity 



1 

t7 



<io-''s- 




FIG. 1: The single electron spin relaxation rate in a GaAs QD 
versus magnetic field (according to Eq. 1351 '). derived for the 
hyperfine-phonon mechanism discussed in the present article. 
The solid line corresponds to the case when the confining 
potential vibrates together with the lattice, while the dashed 
line refers to the case of a fixed confining potential. 



of transverse sound waves, s ~ st ~ 3 x 10"^ m/s^^. The 
typical transverse dimension of a quantum dot is zg ~ 10 
nm and its lateral size is I ~ 100 nm. The dot contains 
about N ~ 10^ unit cells (8 nuclei in each). Hence, we 
can write for the relaxation rate the approximate expres- 
sion: 



— ~ 1 X 10"^^ — 
Ti 1 



UJ 



■f{ki). 



(35) 



For the transition between Zeeman energy sublevels 
the phonon energy is equal to the Zeeman splitting, 
\g*^BB\, which corresponds to 0.025 meV T~^ in GaAs 
where g* = —0.44. At the same time the cyclotron en- 
ergy, hcuc, is as much as 1.76 meV T^^ due to the small 
effective mass in GaAs, m* = O.OGTwe. Consequently, 
the condition lu <ti wq is always satisfied when a perpen- 
dicular magnetic field is applied. For in-plane magnetic 
fields a; < Wo when i? < 10 T in a lateral and B < 100 
T in a vertical QD (the single particle energy spacing in 
a lateral QD is 100 - 300 fJ,eV, and in a vertical QD the 
confining energy of an approximate 2D harmonic poten- 
tial is ^ 4 meV^^''^^). The condition fc^ ^ 1 corresponds 
to a magnetic field of the order of 1 T. Thus, for mag- 
netic fields larger than ~ 1 T the relaxation rate is linear 
in the emitted phonon frequency and as a consequence 
it is linear in the strength of the magnetic field (see Fig. 
([T|) which is plotted for the limit kzo < 1 corresponding 
to _B < 10 T). For the fixed or "independent" confining 
potential we would have f{kl) = 1 and the dependence 
of the rate on the emitted phonon frequency would be 
linear even for small magnetic fields, i.e. when fcZ <C 1 
(apart from the temperature-dependent factor (n^ + 1)). 
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As to the order of magnitude, the mechanism that has 
been considered in this paper gives a rate of the order 
of lO"** s~^ for a temperature of about 1 K and a mag- 
netic field of about 1 T. It therefore appears to be much 
less efficient than the piezoelectric coupling mechanism 
considered in Rcf. 19, where the rate is about 1 s^^ for 
comparable values of temperature, 4 K, and magnetic 
field, 0.5 T. 

For the triplet-singlet transition in a lateral dot we 
find a rate of the order of 10~* s^^ for temperatures up 
to ^ IK {hoj ^ lOO^eV), and in a vertical QD we have 
l/Ti ~ 10""^ s~^ for temperatures up to ^ 10 K and 
hoj ^ 1 meV (we note that in both cases kl ^ 1). The 
latter result should be compared with 1/Ti « 2 x 10^^ 
s~ calculated in Ref. 

Let us now turn to spin relaxation in silicon. Taking 
into account the natural abundance of ^^Si nuclei with 
non-zero spin = 4.68%, their magnetic moment 

jjLj = — 0.56/iAr, the lattice constant a — 5.43 A and the 
electronic density at the position of the nucleus r] ~ 186^ 
we find the effective hyperfine coupling constant to be 
A ~ 5 /xeV. This is far smaller than in GaAs, due to 
the smaller j] and smaller percentage of nuclei with spin. 
Inserting the mass density of Si, p ~ 2.3 x 10'^ kg/m"^ 
and the transverse sound velocity St — 5.4 x lO'^m/s, 
we obtain a prefactor of the order of 10"^" in Eq. (|35|l . 
This corresponds to a very small relaxation rate l/Ti ~ 
10~^ s~^ between Zeeman sublevels in a magnetic field 
of _B = 1 T, for temperatures up to 1 K (g = 2 in Si). 
Note that the spin-orbit related electron spin relaxation 
time in lateral Si QDsi^, for Si:P bound electrons and 
for QDs in SiGe heterostructures^'^ has been predicted 
recently to be on the order of several minutes for the 
same values of magnetic field and temperature. In the 
latter case, however, the relaxation rate may be strongly 
decreased, below the values found for our mechanism, by 
application of uniaxial compressive straini^. 

In general, electron spin relaxation induced by 
electron-nuclei hyperfine interaction in a quantum dot is 
not as efficient as relaxation due to spin-orbit interaction 
for typical values of system parameters^iiSiiiiiSiiSiiiii^. 
In many experiments, the hyperfine related rate would 
be masked by the spin-orbit mechanism (However, see 
Ref. t\2k for ^ case where the hyperfine mechanism 
may dominate). The present experimental data for the 
triplet-singlet transition rate in a GaAs vertical QD is 
Ti « 200/iS at temperatures up to 0.5 K and triplet- 
singlet energy splitting es-t ^ 0.6 meV— . For the re- 
laxation time between Zeeman sublevels in a lateral GaAs 



QD only a lower bound is available: Ti > 50 /xs for an 
in-plane magnetic field B=7.5 T at T=20 m¥M. Both 
of these measurement results were obtained by means of 
transient current spectroscopy^^. 

On the other hand, hyperfine coupling mechanisms 
(such as the one considered in this paper) may be par- 
ticularly relevant as far as effects like dynamic nuclear 
polarization are concerned, where the electron-nuclei hy- 
perfine interaction plays a crucial rol o"^i^^i^^ . Recently, a 
hyperfine nuclear spin relaxation time on the order of 10 
min was measured in a single GaAs QD, at a bath tem- 
perature of 40 mK and a magnetic field up to 0.5 T— . 
However, this experiment dealt with a nonequilibrium 
transport situation, with a resulting spin-flip mechanism 
whose rate turns out to be orders of magnitude larger 
than the one discussed in the present article (see Ref. 

113). 



VI. CONCLUSIONS 

In summary, we have considered a specific mechanism 
for inelastic electron spin relaxation in a QD induced by 
the electron-nuclei hyperfine interaction in combination 
with lattice vibrations. We have found that the inter- 
ference between first and second orders of perturbation 
theory is essential for a correct description of the suppres- 
sion of relaxation at small transition frequencies. The 
relaxation rate has been calculated both for the transi- 
tion between Zeeman sublevels of the orbital ground state 
and for the triplet-singlet transition. We have obtained 
estimates based on these general results and realistic sys- 
tem parameters. The estimates demonstrate that the re- 
laxation rate due to this particular mechanism is very 
small: For the spin relaxation between Zeeman sublevels 
it is much less than the rate calculated earlier in second 
order perturbation theory with an emission of a phonon 
through piezoelectric electron-phonon coupling in a GaAs 
QDiS, and it is less than (or at most comparable with) 
the relaxation rate due to the change in the localized 
electron effective mass induced by the lattice dilation in 
siliconi&. For the triplet-singlet transition, the rate in 
GaAs QDs is still smaller by an order of magnitude than 
the corresponding rate of the piezoelectric mechanism. 

We would hke to thank B. L. Altshulcr, D. V. Averin, 
C. Brudcr, E. V. Sukhorukov and A. V. Chaplik for useful 
discussions. 
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